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ABSTRACT 

We revisit the interlayer tunneling theory of high temperature superconductors and formulate it as a mechanism 
by which the striking systematics of the transition temperature within a given homologous series can be under- 
stood. We pay attention not only to the enhancement of pairing, as was originally suggested, but also to the role 
of competing order parameters that tend to suppress superconductivity, and to the charge imbalance between 
inequivalent outer and inner Cu02 planes in a unit cell. Calculations based on a generalized Ginzburg-Landau 
theory yield results that bear robust and remarkable resemblance to experimental observations. 

Keywords: Superconductivity, high temperature superconductivity, multilayer cuprates, interlayer tunneling, 
competing order, d-density wave 

1. INTRODUCTION 

Prom the point of view of a theorist one of the hardest quantities to calculate is the transition temperature of any 

phase transition, because it is non-universal and depends on many details. Nonetheless, the fascination for room 
temperature superconductivity has been a major driving force behind superconductivity research. Thus one might 
wonder if a theorist has any role whatsoever to play in this affair and one may be justifiably skeptical. A theorist 
often faces two classes of questions: (1) Why docs the substitution of La by Y in a cupratc superconductor raise the 
transition temperature from 30K (LBCO) to 90 K (YBCO)? (2) What separates conventional superconductors 
from high-Tc cuprate superconductors? The answer to the first question may be far more difficult than the 
second. You might also ask why changing D (dcoxyribo) to R (ribo) in DNA (deoxyribonucleic acid) changes 
its properties dramatically? DNA can replicate but RNA cannot. The sensitivity to chemical elements is not 
special to high-Tc superconductors. The Tc changes from 16 K to 9.6 K going from NbN to NbTi. This is 
very typical of emergent phenomena; the properties of a complex system arc more than the sum of its parts. 
Although in common parlance we often like to label a material by its chemical composition, the actual properties 
are vastly diverse even with respect to small changes. Diversity and specificity are common themes in emergent 
matter. The second question, one hopes, is relatively more tractable. There may be a fundamental shift in the 
mechanism — from phonon mediated superconductivity to something else. And this shift a theorist should be able 
to spot, but not without some clues. These clues come not only from detailed electronic, magnetic and structural 
data, but also from hints of patterns. In the past twenty years, much effort has been invested in unravelling 
details, but little in the patterns or the universalities. The aim of this work is to analyze an universal property 
of high-Tc superconductors. 

The phase diagram of high temperature superconductors shows two striking universalities. First, as a function 
of doping, X, Tc is a dome-shaped curve rising at about 5% and dropping to zero at 30%. The second universal 
feature is its dependence on the number of Cu02-layers, n, within a homologous series (see Fig. 1), as shown in 
Fig. 2. The increase in Tc in going from n = 1 to n = 3, about 35%, may seem very modest but in absolute 
numbers it is about 35 K. Considering that the highest transition temperatTire of the conventional superconductor 
NbaCe is about 23 K, this enhancement is spectacular. It is something that is hard to ignore, although most 
theories of high-Tc superconductors tend to ignore this fact and little experimental attention is paid to multilayer 
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Figure 1. Left: a schematic sketch of a homologous series. OP stands for the Cu02-planes next to the apical oxygen 
and IP stands for the inner planes. Note that starting at n = 5 one can in principle distinguish IP from IP*. In actual 
experiments the doping of IP and IP* are very close. The nominal doping level is denoted by x. 

Figure 2. Right: optimal transition temperatures within a homologous series, HgBa2Ca„_iCu„02n+2+(5, as a function 
of n, adapted from Kuzemskaya et al} 

cuprates. There are exceptions, however.^' The inter layer tunneHng theory was motivated by precisely this 
fact.4 

Why is there a superconducting dome in the first place? A moment's reflection will convince one that 
this is quite unusual from the perspective of the BCS mechanism: any arbitrary attractive interaction leads to 
superconductivity. The key to this puzzle, we believe, is competing order. A weakly interacting Fermi liquid 
has very little ordering tendency, simply because most perturbations around it tend to die out. Once we leave 
the regime of validity of a Fermi liquid, we need to deal with many length scales denoted by their characteristic 
coupling constants. It is almost a truism that in a strongly correlated system any symmetry that can he broken 
must he broken. If we ascribe superconductivity to strong correlation effects, a corollary is that a competing 
broken symmetry state may rear its head. Understanding competing order is a key to understanding high-Tc. 

A recent paper by Chakravarty, Kee, and Volker^ traces the second universality to the combination of 
enhancement due to interlayer tunneling and charge imbalance between the layers as n increases. If we denote 
the Cu02-planes closest to the apical oxygens as outer planes (OP) and the others as the inner planes (IP), the 
charge carriers in a nominally optimally doped sample prefer to reside on the OP's, leaving the OP's overdoped 
and the IP's underdoped.®' ^ The aim of the present work is to follow up on these results in greater depth. 

This paper is not a typical conference paper because it reports a substantial amount of original research not 
found elsewhere. In Sec. 2, we set the stage by a discussion of the conventional Josephson couphng energy. 
Sec. 3 is the central section in which an effective Hamiltonian coupling the planes of a high-Tc superconductor 
is developed. This derivation has not been reported previously. Using the results of Sec. 3, we develop the 
Ginzburg-Landau functional in Sec. 4 and apply it to the topic announced in the title. Some of the results 
discussed here have also not been previously reported. In Sec. 5 we present our final thoughts and in the 
Appendix we discuss a technical derivation. 



2. THE FLY IN THE OINTMENT 



In this section, we shall learn an important lesson from the Josephson coupling energy between two identical con- 
ventional BCS superconductors. Assuming a momentum non-conserving tunneling matrix element Tpq between 
two identical s-wave superconductors, the coupling energy Ai? is* 

AE = Ej{l-cos<j)) , (1) 

where at zero temperature Ej is given by 

i?,; = i^|A|(x|^|, ±.=4ne^N\0){\T^^\'), Rq ^ ^ ^ 6.5Kn, (2) 

where A^^(O) is the square of the density of states at the Fermi energy, and (|Tpqp) is a suitable average of 
the tunneling matrix element over the momentum space close to the Fermi energy. The two important points 
are: (1) the lowest energy is zero when the phase difference, 0, between the left and the right superconductors 
is zero, and (2) the coupling energy is non-analytic in the order parameter, tp, or the superconducting gap, A. 
Although the phase-dependent term lowers the energy, the first term cancels this lowering identically. Interlayer 
tunneling cannot enhance superconductivity in a conventional layered superconductor by directly lowering the 
system energy! What it can do is suppress phase fluctuations, which can indirectly enhance Tc. 

There is a nice microscopic interpretation as to why this is the case, as shown in the Feynman diagrams 
m Fig. 3.^ These dia grams can be calculated in a textbook fashion^ using continuum superconducting Green 
functions without electron-electron interaction. The momentum integrations are replaced by energy integrations 
by making use of density of states; the result is Eq. 1. When (f> = 0, the virtual tunneling of the pairs and the 
single electrons both lower the energy of the system. However, because of the superconducting gap, the single 
electron tunneling does not lower the energy of the system as much as it would in the normal state, leading 
to a net positive term that exactly cancels the contribution from pair tunneling. It is possible that this exact 
cancellation is a result of the simplifications mentioned above. In any case, regardless of whether or not the 
cancellation is exact, there is an important lesson here, which is that virtual single electron processes compete 
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Figure 3. Feynman diagrams for calculating the Josephson coupling energy between two identical superconductors on 
the left (L) and the right (R). The diagram (a) corresponds to single particle hops from left to right and back, and the 
diagram (b) corresponds to pair hopping. Both are virtual processes. The normal component of the superconducting 
Green function has a single arrow and the anomalous component has double arrow. The momentum non-conserving 
tunneling vertices are denoted by T. 

with the pair hopping processes. It is this competition that will be crucial in what we shall have to say later. The 

*The derivation of this formula* in the limit that the coherence length is very long compared to the separation between 
the superconductors and close to the transition temperature of the bulk follows from the Ginzburg-Landau theory. 
^We thank D. J. Scalapino for an illuminating discussion of this point. 



second remarkable point is the non-analytic:ity of the resuh. There is no analytic expansion of the ground state 
energy at T = in terms of the order parameter. It is also easy to show that close to the transition temperature 
Ej DC IV'P, and the analyticity is restored. ^'^ Thus, there is no analog of the Lawrence-Doniach^^ model at T = 0, 
while there is one close to Tc. This non-analyticity is a very special feature of a Fermi liquid. 



3. EFFECTIVE HAMILTONIAN 

In this section we shall derive an instantaneous effective Hamiltonian that describes coupled layered cuprate 
superconductors. This will be carried out perturbatively in the tunneling Hamiltonian, Ht, while, in principle, 
arbitrary electron-electron interaction will be allowed. This is really the key: the strong interactions are not 
treated perturbatively. Unlike the Hamiltonian for the conventional Josephson effect between two macroscopic 
superconductors, Ht will be assumed to conserve momentum. There are no reasons for it to be momentum 
non-conserving, unless we wish to describe the effects of impurities. Surprisingly, this perturbative derivation is 
well controlled in a number of circumstances. It is important to stress, however, that for a normal Fermi liquid 
this perturbation theory breaks down completely due to vanishing energy denominators. 

The derivation of this effective low-energy Hamiltonian is important for arriving at the correct Ginzburg- 
Landau functional at T = 0, which we then use to discuss the universality of the transition temperature. In 
principle, it is possible to use this effective Hamiltonian in a more microscopic language, but owing to the 
complexity of the problem we avoid this path. 

3.1. A pedagogic example 

In the previous section we have noted that in the conventional Josephson effect, the lowering of the energy due 
to the uncertainty principle as manifested in both single electron tunneling and pair tunneling completely cancel 
each other when the phase difference, (j) = 0. This is not a generic result, as we illustrate with an analogous 
simple problem. 

Consider the two-site (i = 1, 2) Hubbard model for which the Hamiltonian is 

H2-site = SQ^c\aC-lG " ^ (cI,<TC2,a + h.C.^ + U ^n^^Jl^i, (3) 

where c| ^ creates an electron at site i with spin a, and n^^ is the number operator for an up spin at site i. When 

the on-site Coulomb interaction, U = Q, the problem cannot be treated in perturbation theory in t because of 
degeneracy (vanishing energy denominators), but an exact solution shows that the symmetric state is lower in 
energy than the antisymmetric state by 2t. In contrast, when C//t — > oo and the system is half full (one electron 
per site) , one should include the on-site interaction in the unperturbed Hamiltonian and treat the tunneling term 
as the perturbation. So, to zeroth order, the state of lowest energy corresponds to each site being occupied by 
one electron and is of energy 2eo- The excited states correspond to one site being empty and the other doubly 
occupied. The leading contribution to the ground state energy is second order in t, because the first order matrix 
element vanishes. The 4-fold degeneracy of the unperturbed ground state can in principle be lifted in second 
order degenerate perturbation theory: 



4r / 1 



0). (4) 



Note that we have made use of the completeness relation and the fact that the tunneling Hamiltonian has no 
diagonal matrix element. As far as the low energy, 4-fold degenerate subspace is concerned, we have identified 
the effective Hamiltonian: 

4t^ / 1\ 
= ifo + f Si • S2 - - j , (5) 

where is the spin operator at site i. The terms generated in second order perturbation theory are analogous 
to the terms generated in the Josephson case. In the Heisenberg term an electron makes a virtual transition 
to the neighboring site with opposite spin and the electron on that site replaces the original electron (exchange 



process). In the second term an electron makes a virtual transition and bounces back (direct process). In the 
singlet state, the energy is 

E, = 2eo-—, (6) 

but in the triplet state 

Et = 2eo. (7) 

because of the cancellation of the two processes. There are two important lessons: (1) whether or not the 
effective Hamiltonian lowers the energy depends on the state, and (2) the lowering is non-perturbative in U. In 
this particular case of a finite system, we could have equally well treated [/ as a perturbation. After resummation 
we will arrive at the same answer. For infinitely many degrees of freedom this is far from being guaranteed. 

3.2. Similarity transformation 

Let the total Hamiltonian be denoted by 

H = Hi+H2 + Ht, (8) 

where Hi and H2 are arbitrary interacting Hamiltonians describing the layers and Ht the tunneling Hamiltonian 
coupling the layers. Using a similarity transformation, it is possible to eliminate the tunneling Hamiltonian in 
favor of processes involving two particles, or a particle and a hole, and, of course, an infinite series of such higher 
order processes. A recent development in this direction is infinitesimal similarity transformations, pioneered by 
Wegner^^ and Wilson and Glazek.^^ This technique is yet to be applied to the present problem. Here we shall 
consider only a one-step similarity transformation: 

H,s = e-'He^ = H+[H,S] + ^[[H, S],S] + ---. (9) 
Here S is an anti-Hermitian operator. If we choose 

Ht + [Hi + H2,S]=0, (10) 

then 

iJeff = Hi + H2+^[Ht,S]+--- (11) 

In the representation in which Hi + H2 is diagonal, 

{n\SH = (12) 

provided the energy denominator does not vanish. As mentioned earlier, we consider a momentum conserving 
tunneling Hamiltonian, where 

ifT = E*^(k)(cScgl + h.c.). (13) 

k,<T 

The superscripts refer to the layer index and k = {kx,ky) is the in-plane crystal momentum. For high-Tc 
superconductors, the tunneling matrix elements between adjacent layers are**' ^^: 

t_L(k) = ^ [cos(A;j;a) — cos(A;j;a)]^ : tetragonal (14) 

tj_(k) = ^ [cos(fcxa) — cos(fcya)]^ cos(^fca;a) cos(^fcj/a) : body centered tetragonal (bet) (15) 

For simplicity, we shall consider only the tetragonal case. The bet case requires similar but more careful con- 
sideration of the coupling constants. For two coupled gapless Fermi liquids, a simple analysis shows that a 

momentum conserving tunneling matrix element leads to divergences in the series resulting from the similarity 
transformation because of the degeneracy of the single particle states between the layers. 



In general, if we denote the degenerate subspace of {Hi + H2) by |0) = |0)i |0)2, we get 



I (0| Ek . (k) U^ Jc'^i + h.c.)\n) P 

where |n) refers to an excited state of the uncoupled system. Since the crystal momentum is conserved, the 
excited states must have the quantum number k, which can be used to restrict the intermediate states. So we 
can write 



^ e„,k(iV+l) + e„,k(^-l) 



n^O,k 



where e„,k(-^± 1) are the excitation energies for (A^ ± 1) particles, assumed to vanish near ka = (±7r/2, ±7r/2). 
The tunneling matrix element in Eq. 14 is sharply peaked at ka = (7r,0) and symmetry related points; it is 

vanishingly small around the nodal regions. Moreover, if it vanishes at the nodes faster than the excitation 
energies, the sum can be restricted near the maxima. We denote such a region by the compact set{k*}. Then, 

'^r- 1(01 Ek' ^*i(k')(cw''^c|'^„ +h.c.)|nk)|2 

n^ofe^Ck.} £n,k(Af + l) + £„,k(Ar-l) 

« {0\Hi + H2\0)~-^ E l(0|E*^(k')(cKcL'l. + li.c.)|nk)|2. (18) 

In the second line we have assumed that the excitation energy does not vanish in the set {k*} and can be pulled 
out of the sum by its average; the excitation energy of the (A'' ± 1) systems are the same modulo corrections of 
order 1 /N. We can now extend the summation region back to the entire Brillouin zone and add the terms n = 0, 
because these represent vanishing contributions. Using the completeness of states, we can write the effective 
Hamiltonian as 



{0\H,s\0) ^ {0\Hi + H2\0) - ^{0\ 



E*-(k)(^^& + h.c.) 

k,(7 



|0), (19) 



where D = 2(e,- 



So far we have assumed that the translational symmetry of the ground state is preserved. Now we exploit fur- 
ther robust symmetry arguments independent of the details of the microscopic physics. Consider the interesting 
case where the ground state breaks both translational and 17(1) symmetries but is a spin singlet. In particular, 
we consider the situation where a d-density wave (DDW) coexists but competes with a d-wave superconductor 
(DSC).^^ The important distinctions between DSC and DDW are useful to keep in mind. In the case of DSC, 
the pairing is between particles, and the symmetry of the orbital wave function determines the symmetry of the 
spin wave function. For DDW, the pairing is between a particle and a hole, and the symmetry of the orbital 
wave function cannot enforce any symmetry requirement on the spin wave function; for example, there could be 
a singlet DDW as well as a triplet DDW.^'^ The respective order parameters for the singlet states are 



lb 

(ck.crC-k.o-') = -{coskxtt- COS kya)6cr -a' ■■'DSC, (20) 



(ck+Q,o-'=k,a') = i-{cos k^a- COS kya)da^a' -DBW, (21) 

where Q = (7r/a,7r/a), ip is complex, but ^ is a real scalar. Note the Kroncckcr symbols involving the spin 
indices and the breaking of time reversal symmetry by the DDW order parameter, which is necessarily complex. 
For illustrative purposes, it is useful to plot the the band structure for the coexisting DDW and DSC states, as 
shown in Fig. 4. We have used a normal state energy dispersion et = — 2t(cos kxa + cos kyo) + At' cos kj-a cos kyO. 
Note the energy gaps around the M point and the doubling of the nodes and the spectra because of translational 
symmetry breaking. 



Figure 4. The band structure of the coexisting DDW and DSC state. The energy is measured in units of t and t' = 0.3. 
The magnitudes of the order parameters are \4i\ = 0.4 and \tp\ = 0.2. The normal state chemical potential was chosen to 
be /u = —1. 



Noting that t±{]<.) = t±{]s. + Q), we find that the effective Hamiltonian is 

ri - Tf J. U ^ ST + CV^2/ t(l) t(l) (2) J2) t(l) J2) (2) 

keR,BZ,cr 

^ '-k.cr '^k,cr'^k,<T'^k,cr + '-k+Q,cr'^k,o-'^k+Q,<T''k,<T J + l^^J 

where RBZ stands for the reduced Brillouin zone, which is half the full Brillouin zone. In the simplest approxi- 
mation all the coupling constants in Eq. 22 are equal, but they need not be in a more refined derivation. This 
is obvious because the excitation energies are in principle different. Since we shall use this effective Hamiltonian 
merely to determine the form of the Ginzburg-Landau functional, this is of no consequence. It is now obvious 

that a competing order parameter, such as DDW, can reduce the enhancement resulting from the Josephson pair 
tunneling. This reduction is the analog of the cancellation found in the simplest BCS example discussed above. 

The important point in the derivation of the effective Hamiltonian is the structure of the tunneling Hamil- 
tonian, which is sharply peaked in the region {k*}, away from the nodal region where it rapidly vanishes. The 
non-vanishing excitation energy leads to a controlled effective Hamiltonian via the one-step similarity transfor- 
mation. In the absence of the special structure of the tunneling matrix element (constant in momentum, for 
example), other physical processes must be invoked to suppress the low energy hops. It may be possible to 
recover the same result from orthogonality catastrophe in a non-Fermi liquid^ ^ or from spin-charge separation, 
at least locally.^® Similar arguments have been also presented in the stripe mechanism of spin gap proximity 
effect.19 

For a non-Fermi liquid, an electron creation operator of wavevector k and spin t acting on the ground state 
creates a linear superposition of states that carry the momentum k and spin f. The important point is that 
the act of inserting an electron into the many-body system cannot be renormalized away by defining a single 
quasiparticle or a single superposition of a quasiparticle and a quasihole. If we like, this is the defining property 
of a non-Fermi liquid. This insertion process must generate a superposition of multiparticle states, such that 
the total momentum and spin for each individual term in the superposition are still k and spin |. But now the 
excitation energy is not uniquely related to k. The states of low energy are likely to contribute little to the sum 
because of the orthogonality catastrophe, which generally leads to ground state overlaps of the form 



(7V|ck.|Ar + l) = e-"0)^=f^^^ 



(23) 



*The original interlayer theory* was based on this scenario. 



The quantity cjo is a low energy eutofF, while lUc is a high energy eutoff of the order of the total bandwidth, and 
a is an orthogonality exponent that can depend on electron-electron interactions . Clearly the overlap vanishes 
as uiQ — > 0, that is, in the nodal region. The above expression is simply the 2;-factor which must vanish at the 
Fermi surface in a non-Fermi liquid state. In the superconducting state, the states far above the gap are similar 
to the non-Fermi liquid normal state. However, in the regions of the Brillouin zone where the magnitude of the 
gap is large, the orthogonality catastrophe will be cut off because coq in Eq. (23) will be of the order of the energy 
gap A, and the overlap factor will be (^)". In the regions where the gap is vanishingly small, orthogonality 
catastrophe will act with full force. 



4. THE GINZBURG-LANDAU ENERGY FUNCTIONAL 

In this section we will try to capture the robust features of the phase diagram of multilayer cuprates using a 
suitable Ginzburg-Landau functional. The advantage of Ginzburg-Landau theory is that it allows us to avoid 
the difficult problem of writing down an interacting Hamiltonian appropriate for the cuprates. Instead, we 
can take our cue from the phenomenology of these materials and base our theory on the existence of specific 
broken symmetries. Working at zero temperature simplifies matters further by limiting the number of adjustable 
parameters that appear in the theory. If we assume that the ground state exhibits competing DDW and DSC 
orders, then the mean-field energy functional for an individual layer j is given by 



a' li^jf + X' I + a 0f + A 0| + c; | </.| (24) 



where tpj and (pj are the superconducting and DDW order parameters, respectively, and A is the area of the Cu02 
plane. At T = the physics of these materials is (3-|-l)-dimensional, so we are justified in neglecting fluctuations. 
To account for inter layer tunneling, we must add terms coupling layers ) and j + I; we will assume that longer- 
range coupling between layers is negligible. From Eq. 22 we see that the superconducting order parameters in 
neighboring layers will couple through a term of the form — pc(V'j 0|+i + c.c), where pc is a positive coupling 
constant this is analogous to the pair-tunneling term of Section 2. Ignoring (for the moment) interlayer DDW 
coupling, we arrive at the Ginzburg-Landau functional of Chakravarty et al.^: 

n 

J^=AJ2['^' + A' IVil' + a</.| + A<^| + 5 |V,f - Pc (V'iV'*+i + cc.)] (25) 

This is a zero-temperature energy functional for a single unit cell with n Cu02 planes. Therefore the parameters 
a, a', A, A', g, and Pc do not depend on T, although they can in principle depend on the index j. Previously^ 
values to these parameters were assigned using two criteria: (1) in the single-layer case (n = 1, pc = 0), the 
phase diagram of Eq. 25 should produce a superconducting dome as a function of doping x; and (2) the transition 
temperature as a function of n should describe an approximately bell-shaped curve peaked at n = 3 (see Fig. 2) . 
As demonstrated in Fig. 5, the first criterion can be satisfied by setting A, A', and g equal to constants and 
choosing a and a' as linear functions of x. 

To fulfill the second criterion, we must find a way to include both interlayer tunneling and the charge 
imbalance between the layers. The latter occurs in multilayer samples with n > 3 because the outer layers 
feel a stronger electrostatic potential from the apical oxygens (Fig. 1). As a result, the inner layers of the unit 
cell become underdoped (with respect to the nominal doping x) and the outer layers become overdoped. If 
the nominal doping is optimal, this causes a competing order parameter to appear on the inner layers, which 
ultimately decreases Tc- Knight shift measurements performed by Kotegawa et al.^^ provide values {xj} for 
optimally doped systems {x ~ 0.2) up to n = 5. Using these values, we can incorporate the charge imbalance 
into our functional through the parameters a = a{xj) and a' = a'{xj), whose doping dependence has already 
been fixed by the single-layer case. Then the only adjustable parameter in Eq. 25 is the strength of the Josephson 
coupling, pc, which can be chosen to optimize the shape of Tc vs. n. Chakravarty et al.^ considered three possible 
indicators of T^. ipmax = max\t{jj\, tp^vg = J2i iV-'il/"-! ^^nd V-'rms = E" iV'j P/^^]^^^! where {ipj} are determined 
by a numerical minimization of Eq. 25 with respect to the 2n order parameters. The calculation was performed 
with an open boundary condition, using the downhill simplex method with simulated annealing; the results, with 
Pc = 0.3, are reproduced in Fig. 6. The model clearly captures the peak of Tc sA, n = 3. 
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Figure 5. T = phase diagram of a single-layer cuprate as a function of doping x. \^p\ is the magnitude of the 
superconducting order parameter and <j) is the DDW order parameter. We obtain a superconducting dome by taking 
a{x) = 27(3; - 0.22), a'{x) = m{x - 0.3), A = A' = 1, and 5 = 1.2. 
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Figure 6. Left: T — superconducting order parameters of n-layer cuprates with no interlayer DDW coupling (J = 0). 
Doping, charge imbalance, and Ginzburg-Landau parameters are the same as in Chakravarty et al.^ We plot three possible 
indicators of Tc. Tc oc V'max (+), Tc oc V'avg (x), and Tc oc tp^ms (*). 

Figure 7. Right: modification of the results in Fig. 6 caused by interlayer DDW coupling J — pc, where pc is the strength 
of the Josephson coupling. 



In the following sections we consider two modifications that can affect the shape of the curve for n > 3. First, 
we examine the effects of interlayer DDW coupling, which were neglected in the original model, but should be 
present according to the effective Hamiltonian derived above. Second, we look at what happens when the charge 
imbalance in the five-layer system is allowed to deviate from the numbers reported by Kotegawa et al.^^ 
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Figure 8. Left: correction to tlie results of Fig. 6 obtained with cliarge imbalance R = 1.8 for n = 5. There is no 
interlayer DDW coupling (J = 0). 

Figure 9. Right: same as Fig. 8, but with R = 2.0 for n = 5. When the charge imbalance is this large, ■0avg and ■i/'rms 
return almost to the single-layer values at n = 5, although V'max remains considerably larger. 



4.1. Interlayer DDW coupling 

The effective Hamiltonian derived in Section 3.2 suggests an interlayer DDW coupling of the form J(/)j (/)j+i, with 
J > 0. This is the DDW analog of the .Joscphson term: DDW order parameters in neighboring layers prefer to 
align "antiferromagnetically." (See the Appendix for a more detailed derivation of this result.) We have therefore 
repeated the numerical minimization^ using the modified Ginzburg-Landau functional 

n 

= T+AY,JM3+i (26) 

i=i 

where J is an adjustable (positive) parameter. In Figs. 6 and 7 we compare the results for J = to the results 
for J = (Pc remains fixed at 0.3). When the interlayer DDW couphng is as large as the Josephson coupling, 
the results for n = 4 and n = 5 change quite noticeably, but the peak at n = 3 is robust. 

4.2. Charge imbalance 

The charge imbalance between inner and outer layers occurs only for n > 3, when the layers become crystallo- 

graphically distinct (see Fig. 1). For n = b it is also possible to distinguish between the innermost layer (IP*) and 
its two neighbors (IP); however. Knight shift measurements'^" suggest that the carrier concentration is roughly 
the same on all three. The charge imbalance can therefore be quantified by the ratio R = xq/xj, where xq and 
X] denote the doping of the outer and inner planes, respectively. If the nominal doping x is known, then these 
quantities are also constrained by the relation nx = 2xo + {n — 2)xi. 

The results of the original calculation (Fig. 6) were obtained using values of R and x reported by Kotegawa 
et aZ.^°: R = 1.14 for n = 3, i? = 1.49 for n = 4, and R = 1.64 for n = 5, with x ~ 0.2. The numbers quoted for 
the three- and four-layer materials were averaged over values measured for several Hg- and Cu-based samples. 
By contrast, the value of R given for the five-layer system reflects a single measurement, on Cul245. A more 
recent experiment on the underdoped five-layer compound Hgl245 has yielded R 3.5 for x 0.12.^^ Although 
a direct comparison of i?-values for underdoped and optimally doped samples is probably misleading, there is 
at least reason to believe that the true value for an optimally doped five-layer material could deviate from 1.64. 
Considering this possibility, we have repeated the minimization of Eq. 25 using somewhat larger values of R for 
n = 5. The results are shown in Figs. 8 and 9. 



5. FINAL THOUGHTS 



Wc have focused on the ground state and have pointed out that, despite competing effects, a systematic en- 
hancement of superconductivity can be a universal feature over a limited range of the parameter space in the 
multilayer cuprates. And perhaps a better understanding can provide important directions of research. At fi- 
nite temperatures, thermal fluctuations bring in a set of additional questions that we have not addressed. For 
example, phase fluctuations of the superconducting order parameter are expected to play an important role in 
underdoped cuprates, where the superfluid density is small and the coherence length is short. Indeed, cou- 
pling between the layers can suppress phase fluctuations and correspondingly increase Tc- This has been clearly 
demonstrated in model statistical mechanical systems. It is an interesting question to ask how the interplay 
between fluctuations and competing order will change the shape of the superconducting dome. In fact, a part 
of the downturn of Tc in the underdoped regime can also be due to phase fluctuations rather than a competing 
order parameter. To add more complexity, it may be necessary to step beyond an order parameter theory and 
assess the role of excitations of nodal quasiparticles in the degradation of Tc. The distinguishing feature of the 
mechanism of the superconducting dome discussed here is unique to high-Tc cuprate superconductors, however. 

We are not entirely sanguine that the data shown in Fig. 2 is the final answer, and it would be important to 
explore further the dependence of Tc on the number of layers n. From our theoretical perspective, it is natural 
to suggest that with increasing n, Tc should saturate to the single layer value and the enhancement is only 
an intermediate scale effect. As we have discussed, with increasing n, the outer planes get overdoped, thereby 
decreasing Tc, and the inner planes get underdoped, which also reduces Tc , as long as we believe in the existence 
of a superconducting dome in a single-plane material. This implies that there must be an intermediate layer for 
which the doping is close to the nominal optimal value, which determines Tc . It is difficult to rigorously argue 
that this should be the same as the single-plane Tc, but we suspect that this is the case. In other words, competing 
order ultimately results in an exact cancellation of the enhancement. Of course, other physical processes may 
also arise. It would be extraordinarily interesting to establish the correct Tc versus n behavior from careful 
future experiments. 

Recent experiments suggest that the competing order that is nucleated in the inner layers is antiferromag- 
netic.^^ We are not entirely sure that this the case. Impurity effects can often result in spurious antiferromag- 
netic order. There are further reasons for believing that the aniferromagnetic order parameter may not be 
the key competing order parameter from the "no-mixing of antiferromagnetism and rf-wave superconductivity" 
experiments of Bozovic et al?^ However, /zSR (muon spin rotation) measurements in multilayer cuprates^^ find 
strikingly strong evidence of antiferromagnetism in the inner planes, which appear, however, to be difficult to 
understand. We suspect the positive muons to sit close to the apical oxygens and thus should result in a sig- 
nal dominated by the closest outer planes rather than the inner planes. It is therefore unclear why the /iSR 
signals reflect magnetic moments with magnitudes almost as large as the undoped antiferromagnet. Clearly, it 
is desirable to explore these questions using many different precise experimental tools, as the enhancement of 
the superconducting transition temperature is one of the most important questions in condensed matter physics. 
From the theoretical perspective, however, there is little doubt that a charge imbalance sets in in multilayer 
cuprates, because this effect is the result of robust energetics involving merely electrostatics.^' ^ 

Note added: W. Zwerger has drawn our attention to a paper by R. A. Ferrell, Phys. Rev. B 38, 4984-4985 
(1988), where it is shown that the root of the exact cancellation of single particle tunneling and pair tunneling 
discussed in Sec. 2 is "Anderson's dirty superconductor theorem." Conversely, this cancellation will fail when 
the assumptions of this theorem fail. The simplest suc;h failure is when the gaps of the two superconductors are 
unequal (proximity effect). Two other recent papers on interlayer tunneling are, J. -B. Wu, M. -X Pel, and Q. 
-H. Wang, Phys. Rev. B 71, 172507/1-4 (2005) and M. Mori and S. Maekawa, Phys. Rev. Lett. 94, 137003/1-4 
(2005). We thank the authors for pointing these out to us. 



APPENDIX A. ANTIFERROMAGNETIC DDW COUPLING 



It is sufficient to consider two layers in a pure DDW (rather than coexisiting DDW+DSC) phase. Suppose that 
Hi and H2 both have the form of a mean-field DDW Hamiltonian^*^ on a square lattice, given by 

-ffoDW = X] [(£lk + e2k - M)4<7Ck<7 + (-Elk + £2k - /^)Ck+Q^^Ck+Q,CT + (« Wk4^Ck+Q,<7 + I1.C.) (27) 
keRBZ,cr 

We assume a simple t — t' band structure, with eik = —2t{cos kxa + coskyo) and e2k = 4:t' cos kxU cos ky a. 
The DDW order parameter enters through VFk = {Wo/2){cos kxa — cos kyO), where Wo is proportional to a 
microscopic energy scale times (j). Employing a Bogoliubov transformation, we can diagonalize Eq. 27 to obtain 
the simple two-band Hamiltonian 

Hnuw= [(^k --«)"L"k<. + (i?k -/^)/3L/3k.] (28) 

keRBZ,<7 

The fermion operators ak^ = Cko- + ivu Ck+Q,<7 and /3ko- = sgn( Wk) (^^k Cko- — iuk Ck+Q.o-) destroy quasiparticles 
in the energy bands = e2k ± (^ik + W^)^/^. The complex coherence factors result from broken time-reversal 
symmetry; their real parts are given by 

^ M. + '*'fl-^'-^.. 

Note that Uk depends on the sign of Wk and therefore on the sign of (f> — this is unimportant for an isolated layer, 
but very important for two coupled layers, as we will demonstrate below. Equally important are the symmetries 
of the coherence factors under the transformation k — > k -|- Q: Wk+Q = — sgn(Wk)t;k and t^k+Q = sgn{Wu)uu:, 
which means that the eigenstates created by ak^q ^ and /^k^g ^ are the same as those with momentum k, up 
to a phase factor. In terms of these quasiparticle operators, the zero-temperature DDW ground state |0) can be 
defined for arbitrary filling fraction by ak^|0) cx Q{E^ — fi) and /3k(^|0) oc Q{E^ — /i). 

Now consider the unperturbed ground state manifold |0) = |0)(i) (8) |0)(2) of two identical layers {\(f)i \ = |(/>2|). 
There are two degenerate ground states: |0) = \F), in which (^i = (f>2, and |0) = \AF), in which (^i = —(f>2- Since 
the DDW excitation spectrum is gapped near the (7r,0) points, the arguments of Section 3.2 still apply. The 
effective Hamiltonian for the low energy sector is 

H^s = Hi + H2-^ (30) 

Generally we expect the degeneracy between \F) and \AF) to be lifted by the third term, so that the true 
ground state is some linear combination of the two. What we need to compute, then, are the matrix elements 
of in this subspace. Exploiting the symmetries of the coherence factors, along with the additional symmetry 
t±{ii -h Q) = i_L(k), it is straightforward to show that 

Ht\0) = E {i^(k) sgn [lyf ] - v^^K'^^] 4l'^/3^'j|0) + [1 - 2, 2 - 1]} , (31) 

keRBZ,<7 

provided the system is at or below half-filling. For |0) — |F), we have v^^'' = v'^^ and u^^'' = u'j^\ which obviously 

gives Ht\F) = 0. This cancellation of coherence factors does not occur when |0) = |AF), because then ^k^' = v'^^ 

and Uk^"* = ^Wk^'. Thus we have arrived at a rather unexpected result: the two states are not mixed at any order 
in perturbation theory because Ht\F) vanishes identically. Instead, the state \F) remains at the same energy 
while the state \AF) has its energy lowered by an amount 

AE^^ = -^^ E t^{^rulvle{E+ - f^)e{„- E^) (32) 

keRBZ 

Here we have dropped the superscripts since and are the same in the two layers. The energy scale D is 
set by some combination of the DDW gap Wq and the bandwidth t. 
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